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Critical properties of quantum Hall systems are afltected by the presence of extra edge channels — 
present, in particular, at higher plateau transitions. We study this phenomenon for the case of the 
spin quantum Hall transition. Using supersymmetry we map the corresponding network model to a 
classical loop model, whose boundary critical behavior was recently determined exactly. We verify 
predictions of the exact solution by extensive numerical simulations. 
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More than fifty years after its discovery, Anderson lo- 
calization [l[ remains a vibrant research field. One cen- 
tral research direction is the physics of Anderson transi- 
tions (AT) 0, including metal-insulator and quantum 
Hall (QH) type transitions (that is, betvifcen different 
phases of topological insulators). Apart from electronic 
gases in semiconductor structures, experimental realiza- 
tions include localization of light [3| , cold atoms Jj] , ultra- 
sound [HI , and optically driven atomic systems Theo- 
retically, the field was strongly boosted by the discovery 
of unconventional symmetry classes and a conrplete sym- 
metry classification of disordered systems 0, 0-01 ■ 

Recently it was realized that AT in systems with 
boundaries may exhibit boundary critical behavior dif- 
ferent from, and richer than, the bulk behavior 



The boundary criticality has served to test conformal in- 
variance at two-dimensional AT 13-1^ and as a strong 
constraint on possible theories of the integer QH transi- 
tion [l^. These works often employed the so-called net- 
work models [l7|,|l8| of AT for numerical studies. Within 
a network formulation the richness of boundary critical 
behaviors relates to the possibility of having multiple 
edge channels at the boundary 12 , [l^, [2^ . Physically, 
multiple edge channels occur in an integer QH system 
whenever the filling fraction exceeds one. Some of the 
results below may directly apply to the physics of higher 
QH plateaus and transitions between them. However, 
our description neglects electron interactions at the edge 
which could be relevant in experimental realizations of 
the QH effect 

In this Letter we study boundary critical properties in 
the presence of multiple edge states at the so-called spin 
quantum Hall (SQH) transition j23j . The correspond- 



ing network model [23|, |2J| enjoys a very special status. 



In the bulk, or with reflecting boundaries, the model in 
its minimal formulation (suitable to describe mean con- 



ductances) can be mapped to classical percolation on a 
square lattice 25|, l26[ . This mapping determines exact 
critical properties at the SQH transition 25-29[. In this 
Letter we demonstrate that extra edge channels can be 
straightforwardly included in the mapping. The result- 
ing classical model is not percolation any more, but can 
nonetheless be formulated loop model. 

Both network and loop models (or percolation) are 
lattice regularizations of field-theoretic descriptions of 
AT in terms of sigma-models on symmetric superspaces 
@j Is l [sol . This connection is thoroughly explained in 
Ref. 13 ll ■ and was recently extensively explored by some 
of us |32l |. Through this connection, complete spectra 
of boundary operators were obtained for the conformal 
sigma-models on superspaces cp^+^^-i|^ with a topo- 
logical theta-angle [33[. In the sigma- model approach, 
the number of extra edges is related to the exact value of 
the theta-angle, which affects boundary (but not bulk) 
properties [3J|. The case N ^ M = 1 is directly relevant 
to the SQH effect, and we here apply results of Ref. [s^l 
to obtain exact exponents describing scaling of the mean 
boundary point contact conductances in the presence of 
multiple edges. 

We also report extensive numerical simulations of 
mean conductances in network models on open strips 
with edge channels on both sides. Conformal invariance 
(which has been numerically demonstrated for this tran- 
sition ) relates the exponential decay of the mean con- 
ductance along the strip to dimensions of certain bound- 
ary operators. We extract these dimensions and compare 
them with the predictions of Ref. js^] ■ 

The network model for the SQH effect with extra edge 
channels is shown in Fig. [T] The bulk of width 2L con- 
tains alternating up- and down-going columns of links. 
In addition, m (n) extra columns with the same chirality 
are added at the left (right) edge. These extra links can 



be directed up or down at either edge. We label the four 
possible variants by (£ = zLm, TZ = ±n); positive labels 
(£, TZ) mean the same direction of the edge links and the 
closest bulk link. The links of the network carry doublets 
of complex fluxes (labeled f, i) whose scattering on links 
is described by matrices uniformly distributed over the 
SU(2) group. The scattering at the two types of bulk 
nodes (labeled S = A or B) is described by orthogonal 
matrices diagonal in spin indices: 
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with ts the strength of the quantum tunneling. The SQH 
transition occurs when tA = ts- The boundary nodes 
where the fluxes on the extra edge links scatter, arc de- 
scribed by matrices of the form ([l} but with two inde- 
pendent transmission amplitudes, and tji, one for each 
edge of the system. 

We employ the supersyinmetry (SUSY) method for 



network models [25|, 
the extra edge links 
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with modifications due to 



The row-to-row transfer matrices 
X and Y (formed by multiplying all node transfer ma- 
trices, Ta or Tg for bulk nodes, and for bound- 
ary nodes, in a given row) act in the tensor product of 
bosonic and fermionic Fock spaces defined for each col- 
umn of links. The columns form sites (labeled i) of a 
one-dimensional quantum system whose evolution in the 
vertical (imaginary time t) direction is given by the oper- 
ator U = Yit ^ {^^)^ where Lt is the number of A nodes 
(or B nodes) along the i-direction (see Fig. [1]). With 
periodic boundary conditions in the i-direction, physical 
quantities, including conductance, may be written as cor- 
relation functions, (...) = STr[ ...[/], and the system is 
invariant under a global sl(2|l) SUSY [1^. 

Averaging over disorder independently on each link 
(we denote such averages by overbars) projects the Fock 
space of bosons and fermions onto the fundamental V 
(dual- fundamental V*) 3-dimensional irreducible repre- 
sentation (irrep) of sl(2|l) on up links (down links) j25| . 
The average node transfer matrices T s^i act in the tensor 
products of superspins Vi®V*j^-^ in the bulk and Vi®Vi+i 
or V* ® V*^i in the extra edge regions. Thus, we have 
four types of supersymmetric spin chains 
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All tensor products between neighboring sites of these 
chains decompose into two sl(2|l) irreps. The averages 
Ta b in the bulk and r at the boundaries read 



TA.^ = {l-t'A)I + t\E,, 
Tl,. = (l-ii)/ + <iPM+l> 



A^ B, (6) 
R, (7) 
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FIG. 1. Left: Network model with bulk region of width 
2L = 4, and m = n — 2 extra channels at the left and 
right boundaries. Right: Schematic representation of the de- 
composition of average node transfer matrices. The two top 
rows show bulk nodes, see while the two bottom ones are 
boundary nodes, see ((T)). 



where / is the identity, Ei projects onto the singlet in the 
decomposition of Vi (g) V*^i (and V* (g) T^+i), and Pi^i+i 
is the graded permutation of states on sites i and i + 1. 
The chain ^ with m ^ n corresponds to a sigma-model 
with topological angle Q = (2n + l)7r and to a transition 
between the n-th and the {n + l)-st plateau. 

The decompositions ([61 [7l) have a natural graphical rep- 
resentation shown schematically (without coefficients) in 
Fig. [T] (right). At a bulk node links can be separated 
into disjoint lines in two ways, while at a boundary node 
the lines can avoid each other or cross. Multiplying the 
transfer matrices to calculate the partition function, the 
result is the sum of all contributions of dense closed loops 
of weight one filling the links of the network, weighted by 
factors of either t'g or (1 — ^1) for each node. In the bulk 
the loops are percolation hulls, and the loop model is 
equivalent to bond percolation on a square lattice [25 1. 
The presence of extra edges generalizes this non-trivially, 
since loops can intersect at the boundary. In this situ- 
ation all configurations of the loop model fall into dis- 
joint sectors |37| labeled by the number £ + 2k oi through 
lines extending throughout the system in the i-direction, 
where £ = \C — TZ\ is the minimal possible number of the 
through lines for given C,TZ. 

In the strip geometry our model has external links at 
the top and bottom, where we can insert or extract cur- 
rent. In this geometry the mean conductance g^'^ be- 
tween the top and bottom contacts is the average number 
of through lines going from the source to the drain of the 
current, times 2 for the spin. With the source at the 
bottom, the minimal number of available through lines 



is fc^j'i^ = max(0, C — TZ), whence 



-c.K 



= 2k^;^ + 2 J2Z^ kP{k, Lt/L, Lt/S), 



(8) 



where P{k, Lt/ L, Lt/S.) is the probability (symmetric 
in C and TZ) that exactly 2k "paired" through lines run 
through the system of size 2Lt by 2L + m + n, and ^ is 
the bulk correlation length. At the transition, ^ = (X), 
and for large Lt/L we expect that P{k,LT/L,0) ~ 
g-TT/i ' {k)LT/L _ Conformal invariance at the transition 



3 



allows to identify the exponents h^'^{k) with dimen- 
sions of certain boundary operators. In the literature 
on self-avoiding walks such exponents arc called water- 
melon exponents, and the through lines are called legs. 
For Lt/L 1 the sum in Eq. ([5]) for g is dominated 
by k ~ 1, and we denote j^'''^ the exponent of the first 
subleading correction: 
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(9) 



Critical properties of the above geometric loop model, 
including the exponents h^'''^{k), can be extracted from 
the anisotropic limit of the spin chains ([JHS]), obtained 
by taking all ts I. In this limit the evolution operator 
in one unit of time becomes XY « exp ( — tAtsH) . For 
definiteness we focus on the case the critical {tA = 
ts) Hamiltonian H is then 



2L+m-2 
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The interaction between the supcrspins in the bulk is an- 
tiferromagnetic and uniform at the critical point. Inter- 
actions between boundary spins are ferromagnetic (com- 
pare with d^), and their magnitudes u = {ti^jtA)^, 
V ~ (tR/tA)'^ are kept as arbitrary positive numbers. 
From the diagonalization of H, described in [s^ along 
with full details of the derivation of the results presented 
below, we obtain the critical exponents as scal- 

ing limits of the lowest eigenvalues in a given sector. A 
more detailed knowledge of g^'^ is needed for predicting 
jC,n Q analytically [381, but in most cases we 

can identify its numerical value as h^'^{2). 

The exact expressions for h"^'"{k) derived in [s^l ap- 
pear in Table H There hr,s = [(3r - 2s)'^ - l]/24, 
m = min(m, n) and the parameter is given in terms of 
TO, n and k for m > n (if m < n exchange m and n) by 



6 / \/3 (m + l + k)(n + l-k) 

rk = — arccos \ / 

TT I 2 Y (m + l)(n + l) 



A few remarks are in order. The exponents are 
independent of the boundary couplings w, v, implying 
a boundary renormalization group flow to the stable 
fixed point with infinite boundary couplings, similar to 
Ref. [39| . Moreover, since u,v are positive (ferromag- 
netic), randomness in u, v does not change the exponents. 
If the number of legs is ^ or to -f n -I- 2j with j > 0, the 
exponents equal those of critical percolation with respec- 
tively or 2j hulls through the system (in particular, this 
is so when to = or ?i = 0). However, when the number 
of legs is between £ + 2 and m + n, exponents are highly 
non-trivial and, remarkably, they are irrational. 

The exponents h^'^{k) are symmetric in C and TZ 
by invariancc of the spectrum under left-right reflection. 
Further exponent relations follow from symmetries of the 
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TABLE I. The watermelon exponents, k is the number of 
paired through lines on top of £ unpaired ones and m — 
min(m,n). /i,.,s is the Kac table and is defined as in 
Eq. dni)- 

critical spin chains ([2H5]). In each case V (E) V* (and 
V* (E) V) interacts through Ei, and V(E)V (and V* (E)V*) 
through Pi, and the u,v can be set to one thanks to 
the universality. Then the top-bottom reflection switch- 
ing V and V* induces a mapping between the chains 
(sometimes with different lengths of the bulk region): 
(£,7^) O (-£ - 1, -7^ - 1). This implies: 



(fc) h 



(12) 



In particular, this relates pairwise exponents for the 
chains (0 El), as well as for ^ [4|: /i-"--"(A:) = 
/i™-i'"-i(fc), /i-™'"(/fc) = /i"-i'-"-i(/s). When the to- 
tal number of legs in the model ([5]) equals m + n + 2j 
with j ^ 0, the corresponding exponent /i™^^'"^^(fc) = 
hi,3+2j = {j + l)(2j + l)/3, since in this case we must 
write the total number of legs as (to— l)-|-(n— l)-|-2j-|-2. 

Finally, consider the chain ([3]) . Regarding the leftmost 
site on the right boundary as part of the bulk, and noting 
that we have identical chiralities at the two boundaries 
(so that £ ^ m + n), the {m + n + 2j)-leg exponents are 
simply hi^2+2j = j(2j + l)/3, j = 0, 1, . . . independently 
of TO and n. Results for ^ easily follow from (|T2|l : the 
(to + 71 + 2j)-leg exponents are /ii,2+2j" = j(2j -I- l)/3 in 
both models (Hg]). 

We now present extensive numerical simulations to ver- 
ify our analytical predictions for the exponents h'-''^{l), 
the independence on tL,R, the symmetry relation (|12p . 
and to determine the subleading exponent 7^'^. We nu- 
merically calculate the conductance of critical SQH net- 
works with extra channels in the strip geometry with 
length 2Lt and width 2L + m + n. We parametrize 
the transmission amplitude in Eq. ([T]) as Ia = sin 9 and 
ts = cos 6. In terms of 9 the SQH transition occurs at 
9 ~ 9c = 7r/4. Similarly, we write tL,B. = sin^e, where 
9,, can be arbitrarily tuned or chosen randomly in [0, tt] 
independently for each boundary node. 

The transmission matrix t for the SQH network model 
in the strip geometry is effectively calculated using the 
transfer matrix method [T^j ; the conductance is given by 
the Landauer formula: g = Ti tt'' . Ideally the exponents 
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TABLE II. Numerically obtained exponents h^''^{l),'y^-^ 
and analytical predictions for h^'''^{l), h^'''^{2) for various val- 
ues of {£, TV), oo means that the numerical exponents are ob- 
tained by finite size scaling [i^l. "rand" indicates randomly 
distributed O^- " — " means that numerical estimates of these 
subleading exponents are unreliable. 

h'--->'^[l) and ^-^■■'^ would be obtained by fitting the data 
to Eq. However, finite-size effects hamper such anal- 
ysis in actual simulations unless 9e = Oc or 9e is random. 
Therefore, we apply finite size-scaling analysis to sys- 
tems of various widths ji^]. We only simulate networks 
corresponding to chains ([21 IH E]) with C < TZ, so that 
^mirf — 0. Table |TT] summarizes our numerical results, 
and includes analytical predictions for comparison. 

First we focus on systems ([2]) with £ = 0,72. > 0. As 
shown in Table HH numerically obtained /i°''^(l) agrees 
with analytical results and ^^''^ coincides with h^'^{2) 
for various TZ and 9e- This suggests that the level-one 
descendant of h°^'^{l), /i°''^(l)-(-l=4/3, does not con- 



tribute to the conductance in these systems. In addition, 
/i°'^(l) = 1/3 and 7'^'"'^ = 2 are verified even for random 
^e, confirming our expectation that randomness in the 
boundary couplings is irrelevant. 

Next, we consider cases C^TZ > 0, where /i^'^(l) is 
irrational. Since the L dependence is weak for 6^ = 9c 
and random 9^, we extract several h™'"'{k) for these 9^ 
without finite-size scaling analysis. We confirm that the 
numerically obtained exponents and j™-"- agree 

well with the analytical and /i™'"(2), respec- 

tively. Note that there are no analytically predicted ex- 
ponents with values between /i'"'"(l) and /i"'^"(2). 

Furthermore, we study the cases £, 7?, < 0, whose ex- 
ponents h'-'--^{k) arc related to those of systems with 
£,7^ > by Eq. (fT^ . Comparing results for (— m, —n) 
and (n — 1, m — 1) in Table [III (and assuming the £ -f^ 72, 
symmetry), we sec that Eq. ([T2l) for k = 1,2 is verified 
even if 9e,^9c. Finally, when £ < 0, 72 > in Eq. @ our 
numerics confirm /i^'^(l) = 1 for various cases. Results 
for the subleading exponent ^'-■^'^ are presently inconclu- 
sive because of large numerical errors. 

Our comprehensive numerical simulations confirm the 
exact analytical predictions for the leading exponent 
and, in most cases, allow us to identify the first 
subleading exponent 7-^'^ as h'-'--^{2). 

In conclusion, we have considered the SQH transition 
with extra edge channels. We have mapped the corre- 
sponding network model to a classical loop model, whose 
boundary critical exponents have recently been obtained 
exactly. Using the mapping we obtain exact critical expo- 
nents at the SQH transition from the exponential decay of 
the mean conductance in the strip geometry. Our exten- 
sive numerical simulations confirm the analytical results 
for the boundary exponents. The demonstrated influ- 
ence of extra edge channels on boundary critical behavior 
should be broadly applicable to other QH transitions. 
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SUPPLEMENTAL MATERIAL FOR "EXACT 
EXPONENTS FOR THE SPIN QUANTUM HALL 
TRANSITION IN THE PRESENCE OF 
MULTIPLE EDGE CHANNELS" 

In this supplemental material, we explain the details 
of the numerical finite-size scaling analysis that we use to 
calculate the leading and the first subleading exponents, 
h'~'^{l) and 7'^''^, respectively, from the conductance of 
the spin quantum Hall network model with multiple edge 
channels in the strip geometry. 

Conductance in the strip geometry 

The ensemble averaged conductance g'''^ is predicted 
(see Eq. (9) in the main paper) to decay as 

(1) 

where 2Lt and 2L denote the system's length and 
the bulk width, respectively. Therefore, the exponents 
h'~'^{l) and ^-^-''^ are numerically calculated by fitting 
the averaged conductance .g'^'^ to Eq. ([T]) as a function 
of the aspect ratio Lt/L. In the actual numerical sim- 
ulations, we typically vary Lt/L in the range [2,40] for 
systems for which the predicted exponent h^'^{l) < 1/3, 
and the number of disorder realizations is then taken to 
be 0(10^) - 0(10^). On the other hand, since the con- 
ductance rapidly decays for systems with the predicted 
value h^'^{l) > 1/3, the range of Lt/L for such systems 
is varied between 1 and 4, and the number of samples is 
taken to be 0(10^) ~ O(IO^). As we have mentioned in 
the main paper, in our numerical study we only consider 
the cases of systems for which the minimal number of 
available through lines is k^^^^ ~ 0. 

Figure [T] shows the Lt / L dependence of the averaged 
conductance g'^'^ for several systems with different sets 
of C and TZ as well as different values of the width and 
9e- Although Eq. ([1]) predicts that the conductance de- 
cays exponentially with the universal exponent /i'^'^(l) 
in the long strip geometry {Lt S> i), the averaged con- 
ductance calculated numerically exhibits an exponential 
decay with exponents that depend on the width and 
except for the cases of £ = 0, 7?. = 1 and C = Tl = 1 with 
9p = Oc or randomly distributed Oe- 

To determine the width dependence of the exponents, 
we calculate h^'^{l) and 7'^'^ by fitting the conductance 
for each width to Eq. ([1]) for systems with various sets of 
{C,TZ) and several 6^- Hereafter, in order to avoid con- 
fusion, we introduce the notations h^' (1) and 7^^' for 
the exponents determined in this way, to stress the effects 
of the finite width L on h'~'''^{l) and 7'^''^, respectively. 
We also introduce a fitting parameter fc^' replacing the 
theoretical value = of the conductance of an in- 

finitely long strip. In each case we ensure the goodness of 



TABLE I. A list for the exponents /if ''^(1) and 7f •'^ calcu- 
lated by fitting the conductance of the finite width L to Eq. 
((!]) for various {C,TZ) and 6^. 
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fit to be almost one and confirm that kj^' is almost zero. 
The fitting results are summarized in Table [J Figure [2] 
highlights the width and 6^ dependence of the exponents 
/i^'^(l) and /i^'^(l). As shown in Fig. [2] (a), even in the 
case of £ = and TZ = 1, h'^£^{l) depends slightly on 
the width for all tested 9^, although the exponent seems 
to approach the analytical prediction /i°'^(l) = 1/3 with 
increasing width L. 

In contrast, h^£^{l) is more sensitive to the value of 
9e as shown in Fig. [2] (b). For Oe < Oc, ^l^^I) shows a 
strong width dependence upon approaching the analyti- 
cally predicted h^-^{l). For Oe > 0^ h]^^{l) is close to the 
analytical prediction even for small widths. Remarkably, 
the width dependence is significantly diminished when Oe 
is set as Oe or randomly distributed at each node. This 
diminishing is confirmed for other systems with C,TZ > 0. 
We note that 7^' also shows the width and Og depen- 
dence in some cases, as found from Table IH 

Summarizing this section, we have found that the ex- 
ponents /i^' (1) and 7^' calculated from the conduc- 
tance in finite size systems depend on the system width, 
and differ from the analytical predictions due to finite- 
size effects. 
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FIG. 1. The averaged conductance g^'^ as a function of Lt/L for systems with different width 2L and 9e, as indicated in the 
figure, in the cases of (a) {£,, TZ) — (0, 1), (b) (1, 1), (c) ( — 2, —2), and (d) (—1, 0). Note that the error bars are smaller than the 
symbols. 



Finite size scaling analysis 

Here we explain how we use a finite-size scaling analy- 
sis to calculate h^'^{l) and j^''^ in the thermodynamic 
limit from the conductance of finite systems. The fact 
that hj^' (1) and 7^' depend on the width L implies 
that there are irrelevant fields at the critical point of the 
SQH network model with extra edge channels. Here, we 
generalize Eq. ^ to include contributions of the leading 
irrelevant field as follows. 

i. In the standard finite-size scaling analysis of the An- 
derson transition a Lyapunov exponent describ- 
ing the exponential growth or decay of eigenvalues 
of the transfer matrix in the long quasi-one dimen- 
sional geometry is considered as a quantity described 
by a scaling function. By analogy, we assume that 
/i^'^(l) and 7^'^ in Eq. ^ are also described by 
scaling functions of the width L, such as 



= A \i2L)'/''uoix),i2L)yu,ix) 



1l 



(2L)i/''uoW,(2L)^Ui(x) 



(2) 
(3) 



where uq[x) and 1x1(2;) represent the relevant and ir- 
relevant scaling fields, respectively, and v and y{< 0) 



are their scaling exponents. The variable x repre- 
sents the distance from the critical point. 

Since we focus only on the critical point {x = 0), 
Mo(0) = and ui(0) becomes a constant. Then, the 
scaling functions fh,-^ [0, (2L)*'ui(0)] depend only on 
the irrelevant field ui. We expand them in Taylor 
series up to the order A^'' (iV') as 



TV" 

^[0,(2i:)V(0)]-5]4'(2L)'^^ (4) 

fc=0 

/,[0,(2Lf«i(0)]=^a;j(2Lf^ (5) 



fc=0 



Note that = /i^^(l) = /i^^^(l) and = 7^-^ = 



iii. Since the prefactors Ci,2 in Eq. ^ also depend on 
the width L, we take account of this finite-size effect 
by considering Ci^2 as functions of 1/2L and expand- 
ing them in Taylor series up to the order A^^ {^c)- 
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FIG. 2. The 9^ dependence of (a) and (b) h];^{l) 

for 2L = 30, 50, 70. The horizontal dashed hnes represent 
the analytically predicted values of (a) /i'^'^(l) = 1/3 and (b) 
/i^'^(l) Si 0.037720. The vertical dashed line indicates Sc. The 
symbols at 6^ — represent /i^ 
with randomly distributed Oe. 



"•^1) and fti'^(l) for the system 



According to all these assumptions we write down the 
scaling function for the average conductance as 



AT" 



1=0 



J2Fl\L,LT)+Y.F;'mLT), 

1=0 



exp 



fe=0 



, (6) 



where k^''^, 



7 f^h 



, Ci are treated as fitting param- 



eters. 

In our actual calculations, we set N'^ = N'' 



2 and 



Nq — N2; = 1; there are then 12 fitting parameters. We 
typically computed the averaged conductances for system 
widths 2L = 20, 30, 40, 50, 60, 70 if the predicted /i^''^(l) 
is less than 1; otherwise we computed the averaged con- 
ductances for system widths 2L = 20, 30, 40, 50. By fit- 
ting the calculated conductances to the scaling function 
([S]), we have obtained the exponents h^''^{l) and 7^-''^ 
which are reported in Table II of the main paper. We note 
that the goodness of fit was almost one in all cases. We 
have obtained the irrelevant exponent y « — 1 for systems 
with weak finite-size effects. For systems with strong 
finite-size effects, e.g. {C,TZ) = (1, 1), (—2, —2), (—3, —3) 
with 6*6 = tt/IO, the irrelevant exponent becomes larger: 
y « -0.7 or -0.8. 

To check the validity of our fitting, we have introduced 
the modified conductance defined as 



at; 



-c,n 



g"^"^ - 2k^ 



C.-R 



1=1 



1=0 



/exp 



fc=i 



(7) 



so that only the relevant term in the thermodynamic limit 
is kept in the numerically calculated conductance g^'^. 
We have computed g^^ for each width by using the nu- 
merically calculated fitting parameters, and results for 
several cases are presented in Fig. [3l We thus confirm 
that the modified conductances for all considered widths 
and 9e for each set of £ and 7?. collapse well to a single 
universal exponential function. 
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FIG. 3. The modified conductance g^'^ as a function of Lt /L for systems with different system width 2L and Sg, as indicated 
in the figure, in case of (a) (£, 7^) = (0, 1), (b) (1,1), (d) (-2, -2), and (c) (-1,0). 



